On the pseudogap phase in high-T c superconductors 
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We describe the approach of the superconducting state as a sequence of cross-over phenomena. 
As the temperature is decreased, uncorrelated pairing of the electrons leads to the opening of a 
pseudogap at T F . Upon further lowering the temperature those electron pairs acquire well behaved 
itinerant features at T% , leading to partial Meissner screening and Drude type behavior of the optical 
conductivity. Further decrease of the temperature leads to their condensation and superconductivity 
at T c . The analysis is done on the basis of the Boson- Fermion model in the cross-over regime between 
2D and 3D. 

PACS numbers: 74.72.-h, 74.25.Dw, 74.25.Gz, 74.25.-q 



It is now generally accepted that the onset of supercon- 
ductivity in high temperature superconductors (HTS) is 
controlled by phase fluctuations whilst a finite am- 
plitude of the order parameter is expected to persist up 
to some T*, which can be well above T c . It has also been 
argued j| that the shrinking of the pseudogap phase (the 
interval [T c , T*}) with increased doping should be related 
to an increase of the superfluid density n s rather than 
of the coherence length. In that case, doping of HTS 
should not induce a cross-over between a Bose-Einstein 
condensation (BEC) of preformed pairs and a BCS state 
and hence the opening of the pseudogap should be unre- 
lated to the onset of superconducting fluctuations. These 
expectations seem to be confirmed by specific heat data 
||, transient Meissner screening Q and Andreev spec- 
troscopy H . The picture evolves according to which 
a pseudogap opens up, exclusively driven by amplitude 
fluctuations (dynamical pair-formation, uncorrelated in 
space and not necessarily related to incipient supercon- 
ductivity). At some lower temperature short-range and 
short-time correlations of the electron pairs set in, even- 
tually driving the system into a superconducting state. 

From a theoretical point of view these are questions 
which can be addressed without having to resort to a spe- 
cific microscopic mechanism for electron pairing and can 
be studied on such models as the generalized BCS hamil- 
tonian, the negative U Hubbard model or the Boson- 
Fermion model (BFM). These phenomcnological models 
H which explicitly incorporate pair correlations, capture 
a number of normal state properties of the HTS. 

In a classical BCS type superconductor long range 
phase coherence occurs as soon as short range ampli- 
tude correlations (in form of electron pairing) set in. A 
simple mean field treatment can perfectly describe this 
situation. On the contrary, in HTS amplitude and phase 
correlations are separated and a minimal treatment must 
contain the possibility of allowing for propagating modes 
of electron-pairs on a time scale of the order of the in- 



verse zero temperature gap Q]. We shall here explore 
the intricate relations between the onset of pairing of the 
electrons, their becoming itinerant and ultimately their 
condensation. In order to separate these various features 
in approaching the superconducting state we study them 
in view of a dimensional cross-over such as to capture 
part of the doping induced changes when going from the 
under into the overdoped regime. The dimensional cross- 
over will be controlled by the degree of anisotropy of the 
electron dispersion. The anisotropy in the electric trans- 
port coefficients || lends itself to such a picture. 

This study is based on the BFM, a phenomenological 
model which is particularly suited for that purpose. The 
notion of rather long lived short range pair correlations 
is introduced explicitly into this model by assuming the 
presence of localized tightly bound electron pairs which 
hybridize with pairs of itinerant electrons. We assume 
a system described by effective sites each of which can 
alternatively be occupied by either a bound electron-pair 
or a pair of itinerant electrons, uncorrelated with each 
other. For the present study we shall assume a homoge- 
neous distribution of both the bound electron-pairs and 
the itinerant electrons. It is however feasible that the 
distribution between those constituents occurs in form 
of a phase-separated phase with dynamically fluctuat- 
ing regions of predominantly bound electron-pairs and 
regions with predominantly electrons, leading to "stripe 
phases" . Previous studies of the BFM H , based on con- 
serving diagrammatic approximations let one envisage a 
well defined finite temperature interval for the pseudo- 
gap phase in ID and 2D, while for 3D such a phase JTcj ] 
should be restricted to a very narrow temperature regime 
above T c . This suggests that the nature of the pseudogap 
phase should depend on dimensionality and that the pair- 
fluctuations which control it have, depending on temper- 
ature, different characteristics as far as superconducting 
fluctuations are concerned. We expect that at least two 
energy scales should be involved here, identifiable as two 
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temperatures: T* = T* and T* < T*, 

• Tp corresponding to the opening of the pseudogap 
in the DOS of the electrons due to the onset of 
strong local correlations leading to electron pair- 
ing, uncorrelated in space. 

• Tg corresponding to those local pair states becom- 
ing well defined itinerant excitations due to the on- 
set of temporal and spatial correlations. 

We shall determine the dependence of T F and T B on di- 
mensionality by monitoring the anisotropy of the ratio of 
the electron mass orthogonal to the basal plane to that 
within it, a = mj_/mu in the bare electron dispersion 
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The BFM is then defined by the following Hamiltonian 
H = ^2 Ekct^cuc + E Y^ frq&q + v ^[6 q c k+qi c_ kT + h.c] 

k,<T q k 

with Eq = Ag - 2/j and v denoting the pair-exchange 
coupling constant, c^i denote fermionic operators for 
electrons with spin a and wavevector k and 6 q describe 
tightly bound electron pairs which will be considered as 
simple Bosons. The chemical potential u is common to 
Fermions and Bosons (up to a factor 2 for the Bosons) in 
order to guarantee charge conservation. We examine this 
model within the lowest order self-consistent conserving 
diagrammatic approximation for which the self-energies 
for the Fermions and Bosons are given by: 

v 2 f d 3 k f 
/mE fl (q,w)=-— J — J deImG F (k - q, u - e) 

ImG F (k,e) { tanh (^ji) ~ tanh ^k^r> 



V 1 f d^Q f 

JmE F (k,w) = — J J deImG F (k-q, 



e — lo) 



ImG B (q,e) y tanh (^j^) - cotanh( -'^f"> 



Gp(k, w) = [w — £k — i?eE F (k, ui) — i/mSi?(k, w)]" 1 and 
G_e(q, e) = [w — Eq — ReY, B (q, w) - Hm'S B (q, de- 
note the fully renormalized Fermion and Boson Green's 
functions which have to be determined selfconsistently. 
The above set of equations for the Green's functions is 
solved by a standard iterative procedure. The integra- 
tion over the momenta k|| is carried out by summing 
over a grid of 1025 points in the interval [0, 7r]. The inte- 
gration over frequencies is carried out slightly above the 
real axis i.e., taking to = Reuj + in with n = 0.01 and 
by summation over 2048 real frequencies in the inter- 
val [—2, 2]. Throughout the present work all energies are 
given in units of the bare basal plane bandwidth D = 8iii . 



The pseudogap features which result from the solution 
of the above set of equations depend strongly on dimen- 
sionality. In ID and 2D the Fermionic self-energy is es- 
sentially determined by the term proportional to cotanh 
which leads to the most divergent contribution. For 3D 
the major contribution to the self-energy comes from the 
term proportional to tanh and does not give rise to a 
noticeable pseudogap effect. A pseudogap is only seen 
as we lower the temperature and approach T c where the 
term proportional to cotanh again becomes important. 

The various parameters entering our Hamiltonian are 
determined such as to reproduce certain robust features 
of HTS. Pinning down the number of Bosons n B = (bfbi) 
in this model is not free of ambiguities, given our poor 
understanding of the doping process in these materials. 
For a system such as, for instance, YBCO the number 
of doping induced bound elelctron-pairs could possibly 
be given by half the number of dopant ions 02 _ (1) in 
the chains, thus varying between and 0.5 per effective 
^site. The number of fermions np = ^2 a {cf a Ci a ) is taken 
to be equal to 1 if the boson-fermion exchange coupling 
were absent. We thus obtain a total number of charge 
carriers n tot — np + 2ub which should be contained be- 
tween 1 and 2. As a representative example we choose 
ntot — 1.25. In order to have np < 1 we fix the bosonic 
level as As = 1.1. In order to obtain values for T* of the 
order of a few hundred degrees K we choose v = 0.1. 

The pseudogap manifests itself as a dip which emerges 
in the DOS of the electrons close to the chemical po- 
tential below a certain temperature T F . Tracing the 
value of the DOS at this energy as a function of tem- 
perature permits to identify this characteristic temper- 
ature. The appearence of this pseudogap is linked, as 
we have discussed previously Q , to a breakdown of well 
defined single-electron excitations close to the Fermi sur- 
face. Concomitantly with this trend, two-electron states 
(local electron-pair resonances) emerge and upon lower- 
ing the temperature, acquire itinerant behaviour below a 
second characteristic temperature T B . T B is defined by 
the condition (see Fig.l) that the imaginary part divided 
by the real part of the Boson self-energy is small, i.e., 
7 f(T)/<(T) ee r B (gy,w fl (flf„,T),T)/w B (9,|,T) < 0.1 
for small q vectors such as gii =0.1. Below T B we find 
w B (q\\) = fi 2 q 2 /2m B {T) for q < 0.2) where m B {T) de- 
notes the temperature dependent mass of the itinerant 
two-particle excitations. 

The two-electron states, becoming well defined as the 
temperature decreases, goes hand in hand with a signifi- 
cant increase of the distribution function for those states 
with small wave vectors. It is this feature which leads to 
a finite value of the amplitude of the order parameter in 
the normal state whilst at the same time phase coherence 
is totally absent. The temperature T c at which supercon- 
ductivity sets in is determined by the condition that the 
two-electron excitations condense in a macroscopic quan- 
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turn state. According to the Hugenholtz-Pines theorem, 
it is determined by A B -2^ — i?e£_B(0, 0) = 0. We plot in 
Fig.l these three characteristic temperatures as a func- 
tion of a. Notice the opposite trends of Tp and T c , ap- 
proaching each other as the isotropic limit is approached, 
(a = 10). This is reminiscent of the experimental situ- 
ation in cuprate HTS materials, as we go from the un- 
derdoped towards the optimally doped regime and might 
hence partly be related to such a dimensional cross-over. 
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FIG. 1. Variation of Tp, T B and T c as a function of de- 
creasing anisotropy. In the insert, the effective mass tub (in 
units of m||(T)) of the itinerant in- plane bosonic excitations 
as a function of temperature (solid line) and the ratio of their 
width divided by their real part; for a — 10. 

In order to track such cross-over behavior in the pseu- 
dogap phase, leading from uncorrelated fluctuating elec- 
tron pairs to their itinerant behavior as the temperature 
is decreased, we now show how such features are related 
to precursor Drude behavior of the optical conductivity 
and partial Meissner screening above T c . We encounter 
here physics similar to that of the 2D non-interacting 
charged Bose gas [ pL2| , where a macroscopic number of 
itinerant bosonic states with momenta k < l/£ (£ de- 
noting the coherence length) act together collectively in 
a way similar to the macroscopically occupied condensed 
states with k = in the superconducting phase. The 
optical conductivity <r(q, oS) and the diamagnetic suscep- 
tiblity x(q, u>) are given by the longitudinal, respectively 
transverse part of the linear response to an external vec- 
tor potential A(x, t) depending on space and time. The 
resulting current is given by 

(Ji(x,i)) = -V f d^xdtKijCx-x'^-t^AjCx'^'), 
mf,c J 

Kyfrt) = ie(f)<[?i(x,t), #(0,0)]) - 2m ll 6(x)6(t)6 ij n F (x) 

where n F (x) denotes the density of the Fermions and 
jf (q, t) = J2w 2&iciL_ q / 2 c k+q /2 the Fourier transform of 
their current density. Putting 

Kij{(\,uj) = 6 i: jiu}a(q,,u}) - (q^ - 6 i:j q 2 )x(<l,u) 

permits to extract the optical conductivity and diamag- 
netic susceptibility. The kernel Kij (x, t) is decomposed 



into two contributions. A first one is given by the sim- 
ple bubble for the current autocorrelation function and 
neglecting vertex corrections. This is justified because 
of the strong incoherent contributions of Gf |J- The 
second contribution to this kernel is evaluated in terms 
of the typical Aslamazov-Larkin diagram which takes 
into account the contributions of the itinerant Bosons. 
Due to the intrinsically localized nature of the Bosons 
there are no direct contributions of them to the electri- 
cal current. However since below T B those bosonic two- 
electron excitations become itinerant they will, according 
to this Aslamazov-Larkin mechanism, contribute to sub- 
stantially enhance the conductivity. This is manifest in 
the dc resistivity (see insert of Fig. 2) which clearly shows 
that upon decreasing the anisotropy ratio a, the resis- 
tivity changes qualitatively from an upturn to a down- 
turn as the temperature is lowered. Similarly, the ef- 
fect of the pseudogap and the precursor to a macro- 
scopic quantum state of the two-particle excitations can 
be tracked in the optical conductivity (Fig. 2). By in- 
spection of Fig. 2 we notice that for frequencies above 
~ 0.01 the optical conductivity drops as the tempera- 
ture is decreased below T = 0.03 (indicative of a rem- 
nance of the pseudogap in the DOS of the single-electron 
states), while for frequencies below ~ 0.01 the optical 
conductivity for those temperatures increases (indicative 
of an emerging precursor "Drude" component as we ap- 
proach T c ). The crossover temperature at T =~ 0.01 
is identified as the temperature T B where the electron- 
pairs become itinerant. As the temperature is lowered, 
spectral weight is shifted downwards from the frequency 
regime [0.01 < u> < 0.04] thus enhancing this "Drude" 
component below 0.01. Such a redistribution of spectral 
weight is in fact observed pT| . For frequencies below 
0.001 (~ 240 GHz, taking D = 1 eV) the optical conduc- 
tivity increases substantially as one approaches T c . 

Similarly, the real part of the diamagnetic suscepti- 
bility x'(q, oj) for small wavevectors (q < l/£, £ 
1 / ^/2TOs(T)w,f denoting the coherence length) and fre- 
quencies is negative and small for high temperatures. It 
is given by the usual Landau diamagnetism of the elec- 
trons arising from the first contribution to the kernel . 
For temperatures below Tp (which is 0.01163 for a repre- 
sentative example of a = 10 3 ) we calculate x'(<Z,w), not 
taking into account the phase fluctuations. Close to T c 
(= 0.01127 for a = 10 3 ) x'(<7>0) however is expected 
to diverge as 1/£|- TB , where f? KTB is the Kosterlitz- 
Thouless-Berezinskii coherence length for phase fluctu- 
ations. Our calculated x'il^) increases significantly 
(roughly by a factor 50) over the Landau diamagnetism 
(see Fig. 3). It is given by the Aslamazov-Larkin con- 
tribution of i.e., x'(0,0) ~ ~a,v 4 ( 24mn e * {T y 2 ^f), a 
being a numerical factor of order unity. The expression 
in the bracket corresponds to that known from the 2D 
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Bose gas |12]. It is this contribution which leads to al- 
most complete Meissner screening. A characteristic fre- 
quency scale thus emerges below which phase uncorre- 
cted amplitude fluctuations of the electron-pairs behave 
in some respects as condensed electron-pairs in the super- 
conducting state. For frequencies u> > 10~ 5 we find the 
usual Landau diamagnetism but for uj < 10 -5 we notice 
a saturation of x'(q, w) which decreases as the temper- 
ature increases. It is in this low frequency regime that 
we can extract from the diamagnetic susceptibility an al- 
most complete Meissner screening for strong anisotropies. 
Its physical meaningfulness arises from the fact that for 
such a system there is phase locking over a given finite 
distance which is controlled essentially by the thermal 
wavelength of the itinerant Bosons fl3[| . 



optical conductivity showing Drude behavior. 
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FIG. 3. Real part of the in-plane diamagnetic susceptibility 
for small q and frequency ui — as a function of temperature. 
In the insert, its frequency dependence for different tempera- 
tures, a = 1000. 
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FIG. 2. The in-plane optical conductivity (arbitrary units) 
for a set of different temperatures and a — 10. The inset il- 
lustrates the d.c. resistivity as a function of T, for a = 10 in- 
cluding Aslamazov-Larkin contributions corrections (dashed 
dotted line) and for a — 1000 including (solid line), respec- 
tively excluding them (dashed line) . The arrow on the T axes 
indicates Tt for a — 1000. 



In this Letter we have examined the effect of pure am- 
plitude fluctuations in the pseudogap phase of HTS as 
the dimensionality of the system varies between quasi- 
2D and quasi-3D. We conclude that the doping induced 
cross-over between the under- and the over-doped regime 
might at least partly be related to such a dimensional 
cross-over, clearly indicating a shrinking of the temper- 
ature regime for the pseudogap phase as the optimally 
doped limit is approached. We identified two charac- 
teristic temperatures in the pseudogap phase. The first 
one, Tp, corresponding to the opening of the pseudo- 
gap due to the onset of local uncorrelated electron-pair 
correlations. The second one, Tg, corresponding to the 
onset of itinerant behavior of those electron-pairs, below 
which we expect partial Meissner screening (very similar 
to that expected for the 2D non-interacting charged Bose 
gas, i.e., in the absence of any phase fluctuations) and an 
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